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(Dated: February 1, 2008) 

We reformulate the density matrix renormalization group method (DMRG) in terms of a single 
block, instead of the standard left and right blocks used in the construction of the superblock. 
This version of the DMRG, which we call the puncture renormalization group (PRG), makes 
easy and natural the extension of the DMRG to higher dimensional lattices. To test numerically 
this proposal, we study several quantum mechanical models in one, two and three dimensions. 
In ID the performance of the standard DMRG is much better than its PRG version, however 
for 2D models the PRG is more efficient than the DMRG in a variety of circumstances. In 3D 
the PRG performs also quite well. 

PACS numbers: 75.10.-b, 05.50.+q, 03.65.-w 
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I. INTRODUCTION 

The elucidation of the low energy properties in quan- 
tum two-dimensional strongly correlated systems such as 
Heisenberg, t-J and Hubbard models is known to be one 
of the central issues in modern condensed matter Physics, 
and a variety of approximate methods have been devised 
to treat this type of problems, not always under a con- 
trollable way. It is believed that by solving these issues, 
what is at stake is the solution of the long sought-after 
mechanism underlying high-T c superconductivity in lay- 
ered cuprate compounds, just to mention one of the many 
possible applications. 

After several years of improvements and extensions 
since its formulation M, H, ||, the Density Matrix 
Renormalization Group (DMRG) method has become by 
now a standard numerical tool in the study of strongly 
correlated systems. We shall not dwell upon its defining 
features here, although we will do it in next sections when 
dealing with specific examples. The reader is referred to 
several reviews in the literature, such as ||. ||.||, [J7J. 

The method is specially well-suited when the problem 
under study is purely one-dimensional or quasi-one di- 
mensional (such as ladders [||), although applications 
to small clusters of quantum two-dimensional lattices 
have also been carried out with success [Q, |pT[ , p2| , 

Nevertheless, a formulation of the DMRG method in 
two dimensions on equal footing as it is conceived in the 
one-dimensional case, is still lacking. The performance 
of the method is by far more powerful when the model is 
formulated in a chain rather than in a two-dimensional 
lattice Eq] . There must be deeply rooted reasons in the 
setup of the method to account for this behaviour. 



There have been various attempts to address this kind 
of higher dimensional extensions using several techniques. 
Here we briefly enumerate some of them: i) the original 
"zipping" method for the DMRG in 2D as introduced in 
Jl0| |. This is considered the standard method. It amounts 
to introduce a one-dimensional sublattice inside the orig- 
inal 2D lattice in order to perform the sweeping process 
characteristic of the finite-system method. We shall re- 
view it in Sect. ||; ii) the extensions based on the classical 
Statistical Mechanics applications of the DMRG on 2D 
classical lattices (L7| [[|~8), |R|- These are based on the 
Corner Transfer Method. Recent extensions from 2D to 
3D lattices have been proposed in |2(J, |l[ ; iii) extensions 
i nsp ired in the Matrix Product formulation of the DMRG 
p2[ have been considered i n p3[ ; iv) extensions inspired 
in the superblock method p4j]B. 

Despite all of these proposals, we believe that there is 
still room for improvements regarding this issue. Specifi- 
cally, we want to formulate the DMRG in such a way that 
going from one to two, or higher, dimensions is a natural 
step. We mean by this that all the components enter- 
ing the DMRG need not be modified accordingly with 
the dimensionaliy of the lattice. Quite on the contrary, 
we pursue a formulation in which all the elements of the 
method are independent of the dimensionality. In other 
words, we want a dimension- independent formulation of 
the DMRG method. Moreover, we also envisage the pos- 
sibility of having the errors arising from the truncations 
of the Hilbert space of states kept under control as the 
size of the lattice increases, at least to the extend that 
the errors are also isotropically distributed. 

In this work we have fully developed this desiderata 
for quantum mechanical problems in one, two and three 
dimensions. Quantum mechanical problems have played 
a paramount role in the development of the DMRG M, [| . 



It is well-known by now that it was the analysis of the 
Wilson's RG failure for the simple particle- in-a-box prob- 
lem which led to the invention of the Density Matrix RG. 
Here we want to follow the same approach and we present 
a new formulation of the DMRG in several dimensions 
and check its validity with one-particle quantum prob- 
lems of several types. 

The rationale behind this approach is that if we find a 
new version of the DMRG for 2D quantum lattice prob- 
lems that fails to reproduce the low energy physics of 
the Quantum Mechanics in 2D, then that version of the 
DMRG is doom to failure when applied to more compli- 
cated quantum many-body problems. 

Moreover, as a spinoff of our work, we have devised 
a numerical method to solve the Schrodinger equation in 
several dimensions with an accuracy and efficiency bigger 
than other numerical methods known so far such as exact 
diagonalization techniques. Namely, we can reach lattice 
sizes which are out of reach for exact diagonalization, 
while keeping the same degree of accuracy. 

This paper is organized as follows: in Sect. Owe review 
the DMRG formulation for quantum mechanical prob- 
lems with an arbitrary number of low-energy states kept 
during the truncation process. Then, a new formulation 
using only one block is presented and we compare both 
formulations by stressing their similarities and main dif- 
ferences. In Sect. [II we present an account of numerical 
results to test the performance of the new RG method 
under a variety of circunstances: lattices of several di- 
mensions, discrete Hamiltonians of different types, varia- 
tion in the number of targeted states, etc. Section Sec. IV 
is devoted to conclusions and further prospects. 



II. ONE BLOCK VS. TWO BLOCK 

FORMULATIONS OF DMRG IN QUANTUM 

MECHANICS 

There are several ingredients entering the standard for- 
mulation of the DMRG method 0: the left (B L ) and 
right (Br) blocks which describe the degrees of freedom 
of the system and the environment, respectively; the sin- 
gle sites, one • or two ••, connecting the left and right 
blocks which serve as a kind of probes to test the reac- 
tion of the system degrees of freedom to the coupling to 
the rest of the environment; the universe U, also called 
superblock (SB), which contains the description of the ef- 
fective degrees of freedom of the whole combined system 
of blocks and sites, U = Bl • Br or U = Bl • 'Br, 
at a given step of the RG process; the ground state 
(GS) wave function obtained after diagonalizing the su- 
perblock Hamiltonian of the system H$b and which is 
the so called target state; the density matrix ps of the 
system obtained after projecting the target state down 
onto the Hilbert space of the system states; the sweeping 
process which consists in moving the probe sites •• from 
left to right and viceversa, back and forth through the 



superblock, and it is responsible for achieving the con- 
vergence of the targeted properties of the system, or in 
the RG language, for reaching the fixed point structure 
of the DMRG transformation. 

However, if we want to make a higher dimensional ex- 
tension of the method, we need to abstract its most im- 
portant and relevant pieces from the rest. We can do 
this by setting up the following question: What is the 
distinctive feature which makes DMRG essentially dif- 
ferent from the Wilson RG [£5| and makes it superior 
and so successful? 

The answer is correlation. It is the introduction of corre- 
lations between blocks by means of the superblock con- 
struction and its constant update through the RG itera- 
tions by means of the sweeping process, that allows the 
DMRG to essentially capture the strong quantum fluctu- 
ations present in low dimensional systems. 
Then, a second question arises: What is the role played 
by the density matrix? In a one-dimensional world 
such as a spin chain, one probing site • naturally di- 
vides the chain into two halves according to the scheme 
U = Bl • Br. The correlation between blocks is restored 
by means of the density matrix constructed out of the GS 
of the universe (superblock). What is important here is 
the restoration of correlations between blocks, no matter 
whether it is achieved with the density matrix or with 
another type of construct. 

When we want to extend this construction to higher 
dimensions we inmediately run into some problems. For 
example, in a two-dimensional world such as a quantum 
spin model on the plane, a single probing point • no 
longer divides the world into two halves. If we want to 
recover this splitting we naturally need to substitute the 
point by a line | , in such a way as to divide the universe 
with the following scheme U = Bl\Br. Thus, in ID a 
probing point is needed to do the job while in 2D it is a 
probing line. 

We observe a big difference between ID and 2D, as far 
as the DMRG splitting is concerned: a probing point is 
something manageable for it has a reduced number of 
degrees of freedom while a probing line is not directly 
tractable for it contains a huge set of states. Actually, a 
line poses a non-trivial quantum problem by itself. Thus, 
it seems that we need to stick somehow to the use of 
probing sites instead of lines even if we want to make 
higher dimensional extensions for the DMRG. 

From this preliminary discussion, we arrive at the fol- 
lowing picture for the basic operations in the DMRG 
method: 

Cutting Process: this amounts to splitting the lattice 
into probing sites and blocks so that the superblock is 
made up of these two type of components. 
Sweeping Process: this amounts to moving the prob- 
ing sites throughout the superblock at each step of the 
RG process and updating the content of the blocks in the 
cutting process. 

It is quite apparent that there is a lot of room to imple- 
ment these ingredients. However, we shall present here 



one of such a schemes and we shall do it by comparing 
it with the standard approach of the DMRG in 2D. The 
new scheme will be based on a simpler superblock struc- 
ture made up of only one block and only one probing site, 
namely U = Bu. We shall refer to this site as a punc- 
ture to distinguish it from the usual probing points in the 
DMRG. Thus, we shall also refer to this kind of single- 
block DMRG method as a Punctured Renormalization 
Group (PRG) method. This situation is in sharp con- 
trast with the standard DMRG program which always 
employs two blocks to build the superblock. Later we 
shall see that the number of punctures can be more than 
one. The primary idea behind the PRG version is that 
the block B will give an effective description of all the low 
energy degrees of freedom except for those associated to 
the puncture •. 

Specifically, the problem we want to solve using DMRG 
techniques is the following Schrodinger equation for one 
single particle in several dimensions: 

H^ = E^, (1) 

with the quantum lattice Hamiltonian given by, 



H„ 



2/h 2 + V n 


n = m 


-1/h 2 


|| n — m | 





otherwise 



(2) 



where n, m are vectors of integer components in a square 
lattice of dimension d — 1,2,3; V n is the local poten- 
tial at site n, h is the lattice spacing and . denotes 
the Euclidean norm and we select the nearest-neighbour 
points. This corresponds to a well-known discretization 
procedure of the Schrodinger equation defined in the con- 
tinuum, Hi! = (-A x + V(x))*(x) = E&, with x = ha 
and each of the components ranging asiii = l,2,...,iVj 
\/i = 1 , . . . , d, where Ni is the number of sites in each 
direction of the square lattice. The size Li in each direc- 
tion is Li = xn — x\ — h(Ni — I). Let N be the total 
number lattice sites, i.e., N = Yli=i ^i- The continuum 
limit is recovered as the double limit h — > 0, Ni — > oo 
leaving fixed Li = Ni x h. We shall not be interested 
here in this limit |27| and thus we set h = 1. 

Let us next review the finite-system DMRG algorithm 
which is used to obtain a reduced set of low energy 
states for the Schrodinger Hamiltonian defined in (Q) 
p7| , [ p8[ , pq . For the sake of concreteness, let us as- 
sume that our lattice is two-dimensional. The one- and 
three-dimensional cases will appear as a restriction and 
an extrapolation, respectively, from this 2D case, as far 
as the DMRG implementation is concerned. 

In the following we shall introduce the different compo- 
nents of the PRG method after having reviewed the corre- 
sponding analogous components of the standard DMRG 
program. In this fashion we shall emphasize the main 
similarities and differences between both formulations. 



FIG. 1: Superblock decomposition of a 2D lattice in the 
standard DMRG method. The two central black circles 
represent the probing sites according to the scheme U — 
B L • uB R and S stands for the system block and E for 
the environment block. 
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FIG. 2: Superblock decomposition of a 2D lattice for the 
PRG method. The black circles represent the puntures, 
which in this case are four, and S stands for the system 
block. 



A Superblock Decomposition of the Lattice 

The most common DMRG-decomposition of the lat- 
tice is that of a superblock formed by two blocks and two 
probing sites: U = Bf • •B§_ 1 _ 2 where the subscripts I 
and N — I — 2 denote the number of sites inside each left 
and right blocks such that their sum - plus two - equals 
the total number of lattice sites N. The index I denotes 
the iteration step of the RG process. In Fig. [l] a typi- 
cal superblock decomposition associated to the standard 
DMRG is depicted in two dimensions. 

In the PRG version the superblock decomposition will 
be of the form U = Bi», where Bi stands for a block 
containing N — 1 sites at the step I of the RG process 
and the site will be also denoted as •/. Hence the size of 
the block Bi remains unchanged as we vary I. 

In Fig. S it is shown one of these superblock decompo- 
sitions for a case with 4 punctures. 

A first and crucial difference between the DMRG and the 
PRG is that in the latter there are no distinction between 
left and right blocks which get fused into a single one. 



B Wave-Function Variational Ansatz 



where \Mi) is the normalized vector 



The superblock decomposition of the lattice in turn 
induces a decomposition of the wave-function of the tar- 
geted states associated to the blocks and the sites. To 
be more precise, let us assume that the total number of 
states kept during the truncation/renormalization pro- 
cess is Ne, which includes the GS and Ne — 1 excited 
states. Then, the superblock wave-function \&i(n) at RG- 
step I and lattice point n is split into the following four 
pieces: 



w 



Ni 



:i--p.,)i*) 



(8) 



Ni = y/1- $(.i) 2 



which vanishes at the puncture •;. The PRG ansatz con- 
sists in the following generalization of the state (R) 



|*) =oi|Afj) + o 2 |«i) 



(9) 



*i(n) 



O-Ne + 1 

a N E +2 
^N E 



n = •;+! 

n = »;+2 



E a fi a NE+2+a R? +3 (n) ne B 



N-l-2 



(3) 



where the symbols *i+i and »z+2 denote two points sep- 
arating the blocks Bf and B§_ t _ 2 : {Lf(n)}^ 1 and 
{-R/ Q t _ 3 (n)} Q J 1 are orthonormal basis of states describ- 
ing the degrees of freedom in the blocks Bf and B ] ^_ l _ 
respectively,!. e., 



2- 



( L T\ L l ) = <W - {Rl+3\ R l+s) 



(4) 



and the free unknown coefficients {a a } a= i at this l- 
stage will be determined later on by means of a diagonal- 
ization/truncation process defining the renormalization. 
They are normalized as 



2N E +2 

a|| 2 = £ < 



= 1 



(5) 



Let us now present the PRG version of the SB wave- 
function. To simplify matters we shall first consider the 
case when Ne = 1. Let $(n) be a trial wave function 
describing the GS of the universe U, with energy Eq, = 
($|J?|$). The first step of the method is to puncture the 
universe at a given site, say •/, so that it can be seen as 
the superblock P;», where B\ contains all the sites of U 
but •/, which is labelled by the integer I denoting the RG 
step. To define the puncture operation we shall use the 
projector operator onto the site •;, i.e. 



where a± , ai are variational parameters to be fixed later 
on diagonalizing the SB Hamiltonian. A trivial but im- 
portant observation is that the state |<fr) is in fact a 
particular case of the state |\&) where a,\ — Ni and 
»2 = *£(•;)■ The key idea behind the PRG method is 
the repairing or improvement of |$) at the site •/ by con- 
sidering the dynamics of that site coupled to the rest of 
the universe. 

The puncture method can be inmediately generalized 
to the case with Ne > 1 states. Let us denote by 
{® a } a =i a se t °f Ne orthonormal states, i.e. 



($«|$/3) 



o a ,p 



(10) 



which give a variational approximation to the lowest lying 
states of the Hamiltonian H . We shall also suppose that 
the restriction of H to the subspace expanded by this 
basis is diagonal, i.e. 



(^HW 3 ) =5, 



a,/3 



T7° 



(11) 



where E% are the lowest Ne energies of H at this stage 
of the RG. Projecting out the site •; from the basis 
{3? a } Q Ji, we obtain a set of Ne vectors which are neither 
normalized nor orthogonal. A new set of orthonormal 
vectors {M"} a ^ 1 is obtained by means of the transfor- 
mation 



Ne 



(9=1 



(12) 



where the matrix O a p diagonalizes the scalar product of 
the punctured states (1 — P tl )\<& a ), namely, 



p., |$) = *(.0hi> 



(6) 



where \»i) is the N-component vector with 1 at the posi- 
tion corresponding to the site •;. Using P., the state |<&) 
decomposes as 



|*)=Mi|M,)+4(«0|«j) 



(7) 



Y, °o><*> Pl3 > (So, j, - $ a ' (•i)*? (»i)) = 5 a ,p (13) 



a', /3' 



The analogue of eq. (uti reads 



(14) 



where ' is the inverse of the matrix O. Eq.(|lJ) leads 
to the following PRG ansatz 



N E 



|*) = £ a a \Mn+a NE+1 \. l ) 



(15) 



unknown matrix elements are on the RHS. To be more 
precise, for the standard DMRG method the superblock 
Hamiltonian exhibits the following 4x4 structure asso- 
ciated to the four pieces of the DMRG decomposition of 
the lattice 



which for comparison with the DMRG eq.(0) we write as, 



^N E 



*i(n) 



E^iOaMf(n) neB l 



n 



1N E + 1 

The normalization condition on the a-parameters is: 



(16) 



Hsb(1) 



H l 


H L. l + 1 


Hl»i +2 HhR 


H 'i + iL 


"•t+i'l+i 


"•t+i«i+a ■"•i+iit 


H 'l + 2L 




*l + 2*l + 2 *l + 2R 


Hrl 


H R»i+i 


Hr. 1+2 Hr 



(19) 



Plugging the ansatz (g) into (18) we arrive at the follow- 
ing expression for the various matrix elements in (tL9J) 



a = 



JVe+1 
a=l 



= 1 



(17) 



As in the Ne — 1 case, the wave functions |3> Q ) are a 
particular case of states \& corresponding to the choices 



» 



» 



Oa~' = C~i, a."'' = $"(•;)• The latter vectors are 
orthonormal and expand a hyperplane of dimension Ne 
embedded in the space of SB wave functions which has 
dimension Ne + 1 ■ The PRG method for more than one 
puncture is a straightforward generalization of the one 
puncture case and we do not give the details. 

Both in the DMRG and the PRG we are using the same 
symbol to denote the RG-parameters a, even though they 
come in different number: a (2A r ^+2)-dimensional vector 
for standard DMRG and a (Ne + l)-dimensional vector 
for the PRG. 



C Superblock Hamiltonians 

The superblock Hamiltonians contain the effective de- 
scription of the lattice superblock at a given step of the 
RG process. This is a reduced description of the whole 
Hamiltonian for the original lattice. Thus, the dimen- 
sionality of these SB Hamiltonians is much smaller than 
that of the original Hamiltonian and this makes their di- 
agonalization something manageable. The dimension of 
HgB depends on the number of targeted states and the 
version of the RG method: (2N E + 2) x (2N E + 2) for the 
standard DMRG and (N E + 1) X (N E + 1) for the PRG. 
In order to obtain the matrix elements of the superblock 
Hamiltonian we need to identify the energy associated to 
the original Hamiltonian H in the wave-function ansatzs 
of (||) or (|l^) with the reduced matrix elements of the 
superblock Hamiltonian, namely, we demand 



H^ = (L?\H\L?) 

H. l+1 . l+1 = (•i +1 \H\»i +1 ) 



H 



a, 



P 






H 
H? 



LR 



(L?\H\R? +3 ) 
= (L?\H\'i+i) 



u R'i+i 



(R? +3 \H\' 



i+i/ 



R 

H. 
H?. 



R*l-\-2 



(•i+i\H\»i +2 ) 



(L?\H\^ 
(R?+3\H\' 



1+2) 



(20) 



As it happens, all the matrix elements of the su- 
perblock Hamiltonian have a natural geometrical mean- 
ing in terms of the interactions among the different parts 
of the lattice superblock. For instance, those in the first 
line describe the interactions only within each of the left 
and right blocks; those in the second line represent the 
original interaction at the probing sites *i+i and »;+2i 
those in the third line represent the interaction between 
the left and right blocks and between the two probing 
sites, respectively; and so on and so forth. 
Depending on the type of original interaction, the struc- 
ture of the superblock Hamiltonian will differ. As an 
example, for a problem with a local potential in real 
space, the short-range structure is also transported onto 
the structure of the superblock. Thus, several matrix 
elements are vanishing, like Hlr, H E . l+2 and Hr, 1+1 
fl26| , |27J1 . However, when the problem has long-range 
interactions then all matrix elements are generically non- 
vanishing and the structure becomes more involved 1281 . 

Similarly, to obtain the superblock matrix elements in 
the PRG method we plug ansatz ( |l6| ) into (18) arriving 
at the following expressions 



H SB (l) 



Hr Hb», 

Hr»i a»i 



(21) 



(*i|ff|*i) = (a.\Hs B (l)\&) 



(18) 



Inserting the respective ansatzs into this relation one ob- 
tains an identity between two quadratic forms in the vari- 
able a, such that the input data is on the LHS and the 



H%. t = (M?\H\*) 
H. t = (•i\HM 



(22) 



Similar geometrical considerations apply to the meaning 
of the different matrix elements in terms of interactions 
among block and punctures. 

It is quite apparent that the superblock structure of the 
PRG method is much simpler than the standard DMRG. 
This would imply in principle that the PRG is more eco- 
nomical and advantageous than the DMRG. However, 
we shall see that there is a tradeoff between economy of 
resources and computational time during the sweeping 
process. 

In the DMRG method the entries of the SB Hamilto- 
L9[) are data which are stored as functions of the 



RG step I, and they are updated after every step. How- 
ever in the PRG method it is more convenient to derive 
the SB Hamiltonian from more elementary data. Indeed 
let us consider again the case Ne = 1, where Hb and 
Hb» [ are just two numbers. Using eq.(|8j) we find 



HB = jp ( S * - 2Q{*){»i\H\*) + &{*i){*i\H\*i)) 



(23) 



a entering in the wave- function ansatz (||) and (|l(]), as 
well as the renormalization of the matrix elements of the 
superblock Hamiltonian. 

The starting point is the diagonalization of the su- 
perblock Hamiltonians in order to obtain the wave func- 
tions of the Ne targeted states corresponding to the low- 
est energy eigenvalues. The free parameters a will be 
constructed out of the components of the targeted wave 
functions. 

In the DMRG case, let us denote these parameters as 
the set {a l L ,al [+1 ,al [+2 ,a l R }!*J : 1 , where a l L and a l R are 
TVE-dimensional vectors. 

The truncation of the Hilbert space is performed by 
the projection of the superblock wave- functions onto the 
block formed by one block (left or right) and the near- 
est probing site. Consequently, we have two possibilities 
to perform such truncation, depending on whether we 
project onto the left block or the right block, respec- 
tively, and the explicit form of these truncations reads as 
follows: 
a) B L • •B R -> B L : 



H B . l = iri ((» l \H\$)-*(»i)(»i\H\» l )) 

Hence we can construct the SB Hamiltonian from the 
knowledge of $, its energy E$ and the N-components of 
H\<&), namely (n\H\<f>). Strictely speaking we only need 
(•i\H\4>) in (p3|), but as one wants to move the puncture 
all over the lattice one needs in practice to keep all the 
matrix elements (n\H\</)), Vn. 

Recalling that the state $ is a special choice of the 
general state 9, we derive the important result that the 
lowest energy E of the SB Hamiltonian (|2l]) is lower 
than the energy Eq, of the original ansatz, i.e. 



/ai\ 



a. 

•z + 2 



•i+i 



(26) 



The projected wave functions in the RHS of (pq) must be 
orthonormalized for they will become the new wave func- 
tions of the renormalized left block B' L . Let us denote 
this process as 



= £<? 



a »i+i 



(27) 



Eq < E<$ 



(24) 



which implies that the PRG flow always lowers the energy 
of the ansatz. 

When Ne > 1, using eqs.(|l|) and (p2|) , we obtain the 
entries Hb, Hb» x 



a',/3' 

- * a '(«,)<«i|J3 r |* /9 ') - &\*i){*i\H\$ a ')) (25) 

ffB,rE^ ({• l \H\^)-^(. l ){. l \H\. l )) 



in terms of |$ Q ),i/|$ Q ) and E%. 



D Truncation of Hilbert Space and 
Renormalization 

The next step is to specify the projection of the wave 
functions onto the several blocks and the free parameters 



where O l is the orthonormalization matrix (using a 
Gram-Schdmit method, for instance). If we were only 
targeting the GS wave function, then this would amount 
to a simple normalization of the projected wave function, 
namely, a' L — cll/Ni,, a'. = 0,,/Nl with Nl = \f a \ + a «- 
The new wave functions for the renormalized left block 
L'i +l to be constructed in the next RG-step Z + l are given 
by 



£i*+i(n) = 



E a <a^(n) neB[ 



n 



•l+i 



(28) 



b) B L • •B I 



>B T 



( ^\ 



\*nJ 



(29) 



Likewise, the projected wave functions in the RHS of fl2£ 
must be orthonormalized for they will become the new 
wave functions of the renormalized left block Br: 



The final step of the method is to replace the starting 
wave function $ by the state \& given by eqs. (g) and 
(B4J), while the energy E§ is replaced by the GS energy 
Eq, namely, 






£°"ur 



(30) 



and similarly for the new wave functions for the renor- 
malized right block R'{ +2 '- 



# +2 (n) 



Ea<a-Rz°+3(n) 



n = •z+2 
n e S/+3 



(31) 



Once the truncation process is carried out, we need to 
update the different entries in the superblock Hamilto- 
nian (119). This is an asymmetric procedure, depending 
on whether we are renormalizing from left to right as in 
a), or from right to left as in b). As an example we shall 
give the renormalized left block Hamiltonian in case a). 
In matrix notation it reads, 



H' L (l + 1) = At 



A' 



H L. l + 






A' 



a L 
■,'Ne 



(32) 



(33) 



whereas the right block Hamiltonian H R is plainly taken 
from a previous step corresponding to the appropriate 
length, matching the superblock. These data yield the 
superblock for the next RG-step, namely, B' l I ^ 1 »»B' ] ^_ l _ 3 . 
The case b) of right-to-left renormalization is similar and 
we skip the details. 

Now, let us turn to the PRG method. As it happens, 
things are simpler to formulate in this case: we do not 
need to distinguish between left-to-right nor right-to-left 
cases anymore. There is only one way to truncate. First, 
we diagonalize the superblock Hamiltonian (Ell) and keep 
the Ne lowest energy eigenvectors out of Ne + 1. 
The case Ne = 1 is specially simple since the SB Hamil- 
tonian is a 2 x 2 matrix which can be diagonalized an- 
alytically. The GS energy and wave function are given 
by 



E = - [H B + H„ - ^4 H% mi + (H B - #.,) 2 
H B , t 



ai -D> ° 2 = D(Eo-H mi 



(34) 



D = \ 1 



H B» t 

Eq — H, 



These eqs. are reminiscent of perturbation theory with 
He ( resp. He» 1 ) playing the role of the unperturbed ( 
resp. perturbed ) Hamiltonian. 



|*>-Hfc'> = !*(*!. &)) 

E<$ — ► E^i = Eq 



(35) 



The old ( i.e. $) and new ( i.e.<&') universe wave functions 
satisfy a simple relationship, which can be derived from 
eqs.(||) and (||), namely 



1$') 



N, 



*) 



a 2 - 



Ni 



•i 



(36) 



This eq. means that the wave function $' is a local per- 
turbation of $, accompanied by an global rescaling out- 
side the puncture. For example if the value of $' at the 
puncture is greater than in $, i.e. |$'(»;)| > |$(»;)|, then 
outside the puncture the new wave function will always 
be smaller than in the original one, in order to preserve 
the norm. This interpretation of the PRG method will 
give us some hints to understand the numerical results 
presented in the next section. 

From eq.(]3q) we can also derive the entries of if|$'), 
which will be needed for the next RG step, 



(n\H\&) 



N, 



(n\H\<f>) 



a 2 



Ni 



(n\H\> 



(37) 



Eqs.(p5|), ( pq ) and (p7f) completely define the renormal- 
ized state and we can move to a new puncture, say •;_|_i 
to repeat the process. 

For Ne > 1 we find the lowest Ne states of the SB 
Hamiltonian. Denoting the corresponding wave functions 
as (a 1 , SL) and the energies as E l with i = 1, . . . , Ne, the 
analogue of eqs.(|3q) and (37) read, 



!$'•• 



= £*« 



o aP w 



+ « - Y, 



(n\H\^) = J2 "< O aP {n\H\&) 



(38) 



(39) 



a,P 



+ (a. l -J2*iVc t p$ /3 (»i)){n\H\* l ) 

a, p 

In summary, the PRG method consists in the follow- 
ing steps: puncture the universe at a given site, study 
the block-puncture dynamics and sewing of the puncture 
on the block, getting a new effective description of the 
universe, i.e. 



puncture sow 

U — > B %i — ► 



W 



(40) 



System 



O o^>o 

i i i 

i i i 

o o o 

* * * 

o o o 



o-^o 

I I 



o 
o 



CK> o o^o 



1 1 


o 

1 


[ I 

o o 

o o 


1 



1} 

o 



o^>o CK> o 
Environment 



FIG. 3: Schematic representation of the sweeping process 
in the standard DMRG method. The two black circles 
represent the probing sites which move throughout the 
2D lattice onto the embedded dashed line which serves 
as a guide. 
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E The Sweeping Process 

So far we have defined the RG process at a given step 
for both the standard DMRG an the PRG methods. We 
have then the basic ingredients and we must say how 
the process carries on. Namely, this is the sweeping pro- 
cess which amounts to specify how the original lattice 
is traversed by moving the probing sites and punctures, 
respectively, throughout the lattice. 

In the DMRG case, the sweeping process is graphically 
described in Fig. 0. The probing sites •• move from left to 
right, which enlarges the left block upon renormalization, 
i.e. 



B L 



i m l+i 



r>/L 

B l+1 



(41) 



and similarly from right to left. 

In two dimensions, the sweeping process is acomplished 
by embedding a one dimensional path into the whole lat- 
tice. This line is where the probing sites move along as 
the RG process proceeds. 

In the PRG method the puncture can move through 
the lattice following several patterns. The most conve- 
nient is when each site is traversed 4 times in a cycle or 
sweep: left to right, right to left, up to down and down to 
up. This is graphically depicted in Fig. H Every sweep 
starts at the upper rightmost site and continues with a 
left movement towards the lower leftmost site of the lat- 
tice, following a zigzag line similar to that of the DMRG 
sweeping. At this moment the puncture comes back to 
the original site but following a down-to-up zigzag path. 
So far, each lattice site has been visited twice. To com- 
plete the sweeping cycle, there is another movement, this 
time with a up-to-down zigzag path, towards the lower 
leftmost site. Finally, the puncture comes back to the 
initial site with a zigzag down-to-up path. 
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FIG. 4: Schematic representation of the sweeping process 
in the PRG method. Here the center of the plaquettes 
represent the actual sites of the real 2D lattice. Start- 
ing from the upper rightmost site, the lines indicate the 
successive steps executed to complete a whole cycle or 
sweeping. 

Prior to the sweeping process, there is a warmup pro- 
cess to build up the finite lattice. There is no restriction 
about the initial data of the wave functions, but it is ob- 
vious that a smart guess shall boost the process. In the 
ID and 2D cases we have prepared such a guess using 
a "Kadanoff-like" warmup: the original lattice is coarse- 
grained until the number of degrees of freedom is reduced 
to a given (small) value. Then the effective hamiltonian 
at this scale is exactly diagonalizcd and the resulting 
eigenf unctions are extrapolated to cover the whole origi- 
nal lattice. 

In the 3D case we choose the initial data to be a gaus- 
sian random vector without correlation between different 
cells. Although it takes longer to achieve the desired pre- 



cision, convergence is also completely fulfilled. 



III. RESULTS IN SEVERAL DIMENSIONS 

In this section we present our numerical results. We 
have made an extensive comparative analysis between 
the DMRG and the PRG involving different Hamiltoni- 
ans, lattice's size and dimensionality, number of states 
kept and number of punctures. We consider short-range 
and long-range interactions, in order to sec whether the 
performance of the RG methods depends on this fact. We 
present our results in the following subsections according 
to the dimensionality of the lattices. 

A One-Dimensional Models 

We have solved 3 types of potentials with the two RG 
methods: 

P article-in- a-box Potential: this is the free particle case, 
i.e. V n = 0, Vn. The eigenfunctions are delocalized in 
real space. 

Harmonic Oscillator: this is the standard quadratic po- 
tential V n = n 2 . The eigenfunctions are localized in the 
centre of the lattice. 

2D Delta Potential: this model is a discretization, in mo- 
mentum space, of the two-dimensional delta function po- 
tential Pq| . Due to the rotational symmetry the Hamil- 
tonian is one-dimensional and it is given by 



METHOD 



PBOX a 



OSC" 



flnm — 0n,m9 Q U 



7i-\-7n 



M <n,m<N 



(42) 



where n, m are momentum labels, M and N play the 
role of infrared and ultraviolet cutoffs, respectively, b is 
a control parameter of the discretization which here we 
take as b — y/2 and g is the coupling constant of the delta 
function [§9). 

In Table I we present a summary of our numerical re- 
sults. We have choosen the number of sweeps as the 
basic quantity to establish the comparison between the 
two RG methods. This is because this number only de- 
pends on the RG method employed for each model and 
not on the computer machine. However, to get an idea of 
what this number means we also give the corresponding 
computer time. Moreover, as it happens the time spent 
in a DMRG sweep is different than in a PRG sweep be- 
cause the renormalization operations are also different. 
We have also solved these models using exact diagonal- 
ization techniques. 

Table I shows that the number of sweeps needed by the 
DMRG to achieve the prescribed convergence (10 -10 for 
the Ne/2 lowest energy states) is much smaller than in 
the PRG method. This fact is independent of the model, 
the lattice size, the number of punctures N p and the num- 
ber of targeted states N&- We also observe that increas- 
ing the number of punctures lowers the number of PRG 



Sweeps Time Sweeps Time" 



DELTA C 
Sweeps Time 







0.27 




0.26 




0.02 


G 


2 


0.36 


2 


0.27 


4 


1.5 e 


(2) 


65 


13.92 


37 


7.69 


27 


3.09 


(4) 


48 


4.75 


25 


2.36 


10 


1.25 


(10) 


34 


2.31 


19 


1.25 


4 


0.65 



"Particle-in-a-box potential. 

''Harmonic oscillator. 

C 2D Delta potential. 

d CPU seconds in a Pentium III at 450 MHZ. 

e Here the precision is 10 — 4 . 

TABLE I: DMRG versus PRG numerical results for ID 
systems. For all models Ne = 4 and the number of 
sites is TV = 100 for PBOX and OSC, while N = 38 for 
DELTA. The required precision is 10 -10 for the GS and 
the first excited states. In the PRG method the number 
in parentheses represents the number of punctures N p . 



sweeps, which is nevertheless greater than the DMRG 
ones. On the other hand it is not reasonable to use a 
large number of punctures N p , as compared to the to- 
tal number of sites, for this amounts to an almost exact 
diagonalization of the model. 

The bad convergence of the PRG method in ID is due 
to the "rigidity" of the ansatz, meaning that the up- 
dates of the wave function affects in the same amount 
its left and right handed pieces. The high performance 
of the DMRG can be attributed to the "flexibility" of 
the ansatz, where the left and right handed pieces of the 
ansatz are updated independently. 



B Two-Dimensional Models 

In two dimensions we have solved two models: the free 
particle and the Hydrogen atom whose Coulomb poten- 
tial, in atomic units, takes the following form 



V n = — 



'1Z 



(43) 



y/nJTnl 

with Z the atomic number which is fixed to 1. This 
problem is factorizable and admits an analytical solution 
in the continuum, with energies given by E n — — 4Z 2 /n 2 . 
We adopt a lattice discretization of (|43|) where the origin 
is placed at the centre of a plaquette. 

In table II we display results for a 10 x 10 lattice, which 
show that the PRG is more efficient than the DMRG. 
This fact is independent on the size of the lattice, the 
number of states kept and the model. The aforemen- 
tioned "rigidity" of the PRG ansatz seems to be quite 
appropiate for 2D systems, while the left-right structure 
of the DMRG becomes a handicap. 

The N p = 4 PRG reduces considerably the number of 
sweeps while the consumed (i.e. CPU) time is almost 
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METHOD 



2D-HYDROGEN" 
Sweeps Time c 



2D-PBOX 6 

Sweeps Time c 



DMRG 


9 


45.6 


12 


65.7 


PRG (lxl) 


15 


5.85 


15 


5.85 


PRG (2x2) 


5 


3.23 


5 


3.22 


PRG (3x3) 


3 


3.33 


3 


3.32 



"2D Hydrogen atom 

6 Particle-in-a-box potential in 2D. 

C CPU seconds in a Pentium III at 450 MHZ. 

TABLE II: DMRG versus PRG numerical results for 2D 
systems. For all models Ne — 4 and the number of sites 
is N = 10 x 10 The required precision is 10 -10 for the GS 
and first excited state. In the PRG method the number 
in parentheses represents the number of punctures N p 
(see Fig. |). 



the same as for N p = 9. Both the number of sweeps and 
the CPU time for the PRG are quite insensitive to the 
model. 

In fig. 13 we plot the CPU time tcpu versus the linear 
size of the square lattice Ni — N 2 = L for the H-atom 
( for the free particle we get similar results). From this 
plot we derive the scaling law 



i c pu oc L 



(44) 



Table III collects the numerical values of the exponent 9 
for the Exact, DMRG and PRG methods. They are the 
slopes of the curves in Fig. ||. 



PRG 
DMRG 



FIG. 5: Log-log plot, time icpu versus length size L, 
for the 2D Hydrogen spectrum on a square lattice. Four 
states are targeted in each case, and a relative precision 
of 10~ 10 is required on the ground state and first excited 
state. The PRG uses a 2x2 puncture as in Fig. 2. 



METHOD 


PBOX 


2D-HYDROGEN 


Exact 


6.4 ±0.3 


6.61 ±0.27 


DMRG 


5.5 ±0.1 


5.21 ±0.05 


PRG 


3.92 ±0.1 


3.85 ±0.07 



TABLE III: Values of the exponent 9 in 
with a linear fit of the data plotted in Fig. 



4J) obtained 



Fig. pi shows that the elapsed time is always larger for the 
DMRG than for the PRG. Moreover, the 9 exponent in 
the law (E3) is also larger in the DMRG. This confirms 
that in 2D the PRG method is more efficient than the 
standard DMRG. Table III shows that this conclusion 
remains true for both types of models. For sufficiently 
large lattices the PRG performs better than the exact 
diagonalization methods. 

As an illustration of how the PRG method works dur- 
ing the renormalization process leading to the determina- 
tion of the lowest lying states of the models, we present 
in Figs. 7 and 8 a set of typical snapshots represent- 
ing the four targeted wave functions in a 2D Hydrogen 
model and free particle. We do not plot the three di- 
mensional picture of these wave functions, but instead 
the figure shows their projection onto the x-y plane, so 
that each grey square plaquette is more intense the higher 
the height (in absolute value) of the wave function. The 
punctures of the PRG method (2x2 and lxl) appear in 
these figures as blanck and white plaquettes respectively. 
Making the superposition of all the screenshots of this 
sort, one for each of the RG-steps, one produces a movie 
showing the convergence process of the targeted wave 
functions starting from the warm-up initial states. Dur- 
ing this "time" evolution one sees how the wave functions 
get shaped towards their final exact forms. 



C Three-Dimensional Models 

We believe that the computations presented in one- 
and two-dimensional lattices are enough to make a thor- 
ough comparison of the two RG methods. However, we 
have also carried out calculations in three-dimensional 
lattices. The main purpose of this extension is to test 
the scaling laws obtained in the previous subsection for 
the PRG method. As far as the PRG version is con- 
cerned, there is no need for an extra effort in the compu- 
tation and the sweeping process is a repetition of the one 
employed in two-dimensions, but traversing all the par- 
allel planes making up the whole 3D lattice. However, 
we have not done the similar extension for the DMRG 
method for it was apparent from the previous subsection 
that the higher the dimensionality, the more unnatural 
the standard DMRG sweeping becomes. Thus, we expect 
a worse behaviour in the DMRG method as the dimen- 
sion increases. 



PRG-2 
PRG-3 



1 2.2 2.3 2.4 2.5 2.E 

og(size) 



FIG. 6: Same plot as Fig. H, but for 3D Hydrogen atom. 
PRG is applied with a 2 x 2 and a 3 x 3 punctures. Notice 
the crossing between PRG lines and the exact ones. 

We have chosen the three-dimensional Hydrogen atom 
as the representative model and thus, we have discretized 
the Coulomb potential 



V n = - 



2Z 



\Jn\ + n 



(45) 



with the proviso that the origin is placed outside the 
lattice, at the center of a cube, to avoid the singularity 
of the potential. 



METHOD 



3D-HYDROGEN 



Exact 
PRG 2x2 
PRG 3x3 



9.51 ±0.6 
6.6 ±0.4 
5.6 ±0.3 



TABLE IV: Values of the exponent 9 in 
with a linear fit of the data plotted in Fig. 



4J) obtained 



Fig. g and table IV show that the 9 exponents are 
larger than those in 2D. However, the 3D exponent for 
PRG is again lower than the corresponding one for the 
exact diagonalization method. In fact, the PRG expo- 
nent is not much larger than the DMRG exponent for 
2D lattices. This fact makes us to believe that the PRG 
version is a well-behaved procedure in dimensions than 
one. Moreover, the fact that increasing the number of 
punctures reduces the number of sweeps is also reflected 
in this table by the smaller values of the 9 exponent. 
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motivation to undertake this study is the construction of 
a version of the DMRG which is better suited for higher 
dimensional lattices than the standard version. 
We have stressed the role played by the sweeping pro- 
cess in the DMRG method and we have singled it out as 
one of the most relevant features of the finite-system for- 
mulation of the DMRG. It is this sweeping process what 
becomes one of the defining components of the PRG ver- 
sion of the DMRG and it is responsible for achieving the 
convergence of the lowest energy properties to a certain 
prescribed precision. 

In ID lattices the standard DMRG method outperforms 
the new PRG version. This is natural for we know that 
the DMRG is somehow optimal when dealing with one- 
dimensional models. 

However, in 2D lattices the PRG formulation is more nat- 
ural and well adapted to this type of lattices, unlike the 
DMRG which needs to split the lattice into left and right 
blocks. We have also tested numerically that the PRG 
gives a better performance than the DMRG for several 
types of models. 

We have checked that the PRG method also perfoms well 
for 3D lattices. As a byproduct of this work, the PRG 
version can be considered as a new numerical method 
for solving the Schrodinger equation in 2D and 3D with 
a better efficiency than the exact diagonalization tech- 
niques, as can be seen from Tables III and IV. 

Altogether we find these results quite promising, but of 
course the crucial issue is wether one can generalize the 
PRG to interacting many body systems. The technical 
point is to define a "puncture operation" of many body 
wave functions, which must isolate the "local" states as- 
sociated with the puncture from the "global" states as- 
sociated with the block. As we have shown in quantum 
mechanics this can be achieved by a set of local projection 
operators which "pick up" the value of the wave function 
at any given site ( see eq. |6j). In the many body case one 
needs projection operators for all possible local states of 
the puncture. The formalism that is required to gener- 
alize the PRG to many body systems is reminiscent to 
the Matrix Product (MP) method {30|. This fact may 
not be that surprising since after all the variational state 
underlying the DMRG is an inhomogenous MP ansatz. 
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